The 21cm Signature of a Cosmic String Loop 
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Cosmic string loops lead to nonlinear baryon overdensities at early times, even before the time 
which in the standard LCDM model corresponds to the time of reionization. These overdense 
structures lead to signals in 21cm redshift surveys at large redshifts. In this paper, we calculate 
the amplitude and shape of the string loop-induced 21cm brightness temperature. We find that 
a string loop leads to a roughly elliptical region in redshift space with extra 21cm emission. The 
excess brightness temperature for strings with a tension close to the current upper bound can be as 
high as 1°K for string loops generated at early cosmological times (times comparable to the time of 
equal matter and radiation) and observed at a redshift of z + 1 = 30. The angular extent of these 
predicted "bright spots" is about 0.1° for a value of the string tension equal to the current upper 
bound. These signals should be detectable in upcoming high redshift 21cm surveys. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

21cm surveys promise to become an excellent window 
to observationally probe the "dark ages", the epoch in 
the evolution of the early universe before the main burst 
of star formation sets in which is the time of reionization 
in the standard LDCM cosmology. The reason is that 
the intensity of the redshifted 21cm radiation which we 
receive tracks the distribution of baryons and does not 
require the baryons to have formed stars. 

From an experimental point of view the situation is 
promising. The LOFAR telescope array has been com- 
missioned and has the capability of providing 21cm maps 
between redshift 5 and 100 at excellent angular and red- 
shift resolution [l|. LOFAR can in turn be viewed as 
a prototype for the "Square Kilometer Array" project 
which will have even higher resolution 0, and will also 
be able to probe up to redshifts of 30. 

From the point of view of cosmology theory it is an 
interesting challenge to see whether the current standard 
paradigm of early universe cosmology, the LCDM model, 
will be consistent with the data. In this scenario, the pri- 
mordial fluctuations form a Gaussian random field with 
a scale-invariant spectrum with a slight red tilt (more 
structure on larger scales). In this context, nonlinearitics 
form rapidly on an extended range of scales at late times, 
comparable and later to the time of reionization. Ob- 
serving large amplitude primordial fluctuations at high 
redshifts would pose a serious problem for the standard 
cosmological scenario. 

In contrast, if our theory of matter is of a type which 
admits topological defects such as cosmic strings, these 
defects will lead to nonlinearities at very high redshifts. 
It was recently realized that cosmic string wakes lead to a 
signatures in 21cm redshift maps with a large amplitude 
and with a specific geometrical shape 0). Cosmic string 
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wakes produced at the time t eq of equal matter and radi- 
ation, with string tension close to the current upper limit 
and seen at redshift z + 1 — 30 lead to wedges in 21cm 
redshift maps with an absorption signal characterized by 
a brightness temperature 8Tb ~ lOOmK. The wedges are 
extended in both angular directions (a fraction of one de- 
gree) and thin in redshift direction [8vjv ~ 5 x 10 -6 ). 
Such a signal should easily be detectable to the SKA ex- 
periment. 

It has been realized that many particle physics models 
beyond the "Standard Model" give rise to cosmic strings 
(see [U-[6| for reviews). In models which admit stable 
cosmic string solutions, a network of such strings will 
inevitably form during a symmetry breaking phase tran- 
sition in the early universe. The network will then ap- 
proach a "scaling solution" in which the statistical prop- 
erties of the string network are independent of time if 
all lengths are scaled to the Hubble radius t (where t is 
physical time), and will persist to the current time. Since 
strings carry energy density, their gravitational interac- 
tions will induce cosmological effects. The cosmological 
effects of a string increase in amplitude as the string ten- 
sion (which equals the mass per unit length) [i increases. 
The string tension, in turn, is related to the energy scale 
r\ of the particle physics which leads to string formation 
via ji ~ rj 2 . Hence, searching for cosmological signatures 
of strings is a way to probe physics beyond the "Stan- 
dard Model" of particle physics at the high energy scale. 
It is complementary to efforts to look for physics beyond 
the Standard Model in terrestrial accelerators such as 
the LHC, experiments which probe the low energy end 
of beyond the Standard Model physics. It is possible that 
cosmic strings formed at a very high energy scale are red- 
shifted outside our current Hubble radius by a period of 
cosmological inflation [tJ which occurred after string for- 
mation. However, it has recently been realized that in 
fact cosmic strings are formed at the end of the phase 
of inflation in many models of inflation based both on 
string theory [f| and on more conventional supergravity 
Q. Hence, there is good motivation to search for signa- 
tures of cosmic strings even in the context of cosmological 
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evolution with an inflationary phase There are already 
cosmological limits from the angular power spectrum of 
cosmic microwave background (CMB) anisotropy exper- 
iments which set an upper bound [l2[ (see also [HI for 
earlier work) on the dimensionlcss value G/k of jj, (where 
G is Newton's gravitational constant) which is 2 

G[i < 1.5 x 1(T 7 . (1) 

The scaling distribution of cosmic strings consist of a 
network of infinite strings with mean curvature radius 
comparable to the Hubble radius, plus a distribution of 
loops with radii R < t. A long string segment gives rise 
to a cosmic string wake whose 21cm signals were investi- 
gated in previous work Q. In this paper, we discuss the 
21cm signature of a cosmic string loop. Cosmic string 
loops accrete surrounding matter in a roughly spherical 
pattern [l6[ and lead to overdensities of matter (both 
dark matter and baryons) which form when the string 
loop is produced at very early time and then grow in size 
as time progresses. Thus, string loops will lead to large 
21cm redshift signals at early times (long before the time 
of reionization in the standard LCDM scenario). The 
distribution of these signals is highly non-Gaussian since 
the string loop distribution is not Gaussian. 

The outline of this paper is as follows: In the next sec- 
tion we discuss some preliminaries both concerning the 
distribution of cosmic strings and concerning 21cm cos- 
mology. Section 3 summarizes our computation of the 
brightness temperature pattern induced by an individ- 
ual cosmic string loop. We end with a discussion of our 
results. 



II. PRELIMINARIES 

In this paper we will calculate the 21cm signature of 
a single cosmic string loop. The loop distribution is de- 
termined by the cosmic string scaling solution, according 
to which at all times t (sufficiently long after the phase 
transition which gives rise to the strings) there will be 
a network of infinite strings with mean curvature radius 
comparable to the Hubble radius t. Since the strings 
are relativistic objects, string segments will typically be 
moving with relativistic velocity. Hence, there will be 



1 Cosmic strings may also arise in the form of cosmic superstrings 

10] in alternatives to inflation such as "String Gas Cosmology" 

2 It is important to point out that this limit is sensitive to param- 
eters such as the number of strings per Hubble volume which is 
currently not accurately known. An advantage of looking for di- 
rect signatures of strings in position space (rather than in terms 
of a power spectrum) is that one can obtain bounds on Gfj, which 
are insensitive to the abovementioned paramtcrs. Preliminary 
studies of edge detection algorithms |14j which search for the 
signatures of cosmic strings in position space in CMB anisotropy 
maps have provided evidence |15| that the limits can be improved 
by one order of magnitude. 



frequent intersections of strings. As can be argued an- 
alytically (sec [4-6] for reviews), the scaling solution of 
the infinite strings is maintained by the process of in- 
teractions between segments of infinite strings producing 
loops. 

If the infinite strings have negligible substructure, then 
the loops which they produce at time t will be generated 
with a typical radius of 

Rform(t) = at, (2) 

where a < 1 is a numerical constant (whose value must 
be determined numerically and is at the present time 
still quite uncertain). Once formed, loops will oscillate, 
emit gravitational radiation and gradually shrink in size. 
Numerical simulations of the dynamics of cosmic string 
networks [l7| have confirmed that the loop distribution 
approaches a scaling solution as well, but the detailed 
structure of the scaling solution is still poorly understood 
(see e.g. [l8| for some recent progress). Looking for po- 
sition space signatures of individual string loops has the 
advantage over approaches based on computing correla- 
tion functions that the analysis is much less dependent 
on the unknown parameters entering into the description 
of the string loop scaling solution. 

We will be working in terms of the "one scale model" 
for the string loop distribution (see e.g. [l9j]) which is 
based on assuming that all string loops are born with 
the same radius and that their number density redshifts 
as the universe expands with the radius very slowly de- 
creasing due to gravitational radiation. The result for 
the number density n(R, t) of strings per unit volume for 
unit R is hence 

n(R,t) = vR~ 2 r 2 R>at eq 

n(R,t) = vR- 5/2 r 2 tl/ q 2 -yGfit < R< at eq (3) 

n(R,t) = const R<jG[i, 

where 7 is a constant determined by the strength of grav- 
itational radiation from a string loop and which is in the 
range 7 ~ 10 (with substantial uncertainty). The cos- 
mic string loops which dominate the mass in strings (and 
hence which dominate the baryonic mass accreted onto 
string loops) are hence those with R ~ jG/i. 

Since cosmic strings carry energy, they will have grav- 
itational effects which in turn lead to cosmological sig- 
nals. Most of the signatures of infinite string segments 
originate from the fact that space perpendicular to a long 
straight string segment is conical, i.e. like 1Z 2 with a miss- 
ing wedge, the "missing angle" being 87rG/i [2C|. This 
leads to line discontinuities in CMB anisotropy maps 
pH ] . The cosmological fluctuations induced by a network 
of cosmic strings are "active" rather than "passive" (as 
they are for example in the case of inflationary cosmol- 
ogy). Thus, in spite of the fact that the power spectrum 
of fluctuations induced by strings is scale-invariant (see 
e.g. [lj|), there are no characteristic oscillations in the 
angular power spectrum of CMB anisotropics. The fact 
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that such oscillations have been observed with high sig- 
nificance allows us to put the constraint on the cosmic 
string tension mentioned in ([T]). On the other hand, re- 
cent work [lj| has shown that limits which are up to one 
order of magnitude stronger may be found by analyzing 
CMB data directly in position space making use of edge 
detection algorithms |14| . 

Long string segments moving through the primodial 
dark matter and baryonic gas lead to wedge-shaped over- 
densities in their wake. These are called "wakes" J22j. A 
string at time t will produce a wake which is extended 
(typical length scale t) in the plane of the string world 
sheet, and has a mean thickness of AnGfit. The initial 
density in the wake is twice the background density. The 
wakes will grow in thickness by gravitational accretion, 
and will lead to signatures in CMB polarization [23j and 
in 21cm redshift maps (see also Q). 

The signature of cosmic string wakes in 21cm redshift 
maps is a consequence of the overdensity of baryons inside 
the wake. This leads to extra emission or absorption 
of the 21cm radiation from the direction of the wake 3 . 
The signal of a specific string wake produced at time 
ti is seen at a redshift z(t) corresponding to the time 
t when our past light cone intersects the wake. It is 
extended in both angular directions (with a length scale 
corresponding to the comoving size of the horizon at time 
U), and it is thin in redshift direction direction, with the 
characteristic thickness being Q 



Su 



Hw 



(4) 



where H is the expansion rate of space and w is the wake 
width (taking into account its increase by gravitational 
clustering between the time ti and t), both evaluated 
at the redshift z when the past light cone intersects the 
wake. 

21cm surveys provide a window to explore the dark 
ages of the universe since they are sensitive to the dis- 
tribution of baryons in the universe rather than stars. 
Most of the baryonic matter between the time of recom- 
bination and that of reionization is in the form of neutral 
hydrogen. Neutral hydrogen has a 21cm hypcrfine tran- 
sition line. If we look back at the radiation from the 
direction of a baryon cloud, we will see extra emission 
or absorption in this 21cm line. If the gas temperature 
Tk is higher than that of the microwave photons at the 
redshift z when our past light cone intersects the cloud, 
the signal will be in emission, if Tk < T~ t it will be in 
absorption 4 . 

A significant advantage of cosmology with the 21cm 
line compared to CMB cosmology is that 21cm surveys 



3 It is an emission signal if the gas temperature Tk in the wake is 
higher than the temperature T 7 of the CMB, and is an absorption 
signal if T K < T 7 . 

4 For a detailed review article on the cosmology with the 21cm line 
see El. 



provide three dimensional maps of structure in the uni- 
verse. We map out in two directions in the sky plus 
in redshift direction. The effects of gas clouds at dif- 
ferent redshifts which our past light cone intersects will 
be seen at different frequencies, namely the 21cm fre- 
quency shifted by the respective redshifts of intersection. 
Applied to cosmic string models, whereas the effects of 
different strings which our past light cone intersects are 
all projected into the same two-dimensional maps, the 
21cm signals will remain well separated if the strings are 
at different redshifts. 

In previous work 0] the 21cm signal of an individual 
string wake was studied, and in a followup paper (2(| the 
two dimensional power spectrum at fixed redshift was 
computed for a scaling distribution of string wakes. In 
this paper, we will study the signature of a cosmic string 
loop. 

In the case of cosmic string wakes, we [|| found that the 
signal in 21cm redshift surveys is a wedge in the three- 
dimensional redshift maps which is extended in the two 
angular directions, the angular extent given by the co- 
moving scale of the Hubble radius at the time U that 
the string wake was formed, and narrow in redshift di- 
rection, the relative thickness being proportional to G/i 
multiplied by the linear growth factor z(ti)/z(t) between 
the redshift z{ti) corresponding to the formation time 
and the redshift z(t) at which the wake is being seen 
(at which our past light cone intersects it). For redshifts 
z(t) larger than that of reionization and for values of Gfi 
comparable to or smaller than the current observational 
upper bound (TTJ, the 21cm signal is an absorption sig- 
nal with a brightness temperature amplitude of up to 
400mK, the amplitude being to a first approximation in- 
dependent of Gfi. 

Here, we find that the amplitude of the cosmic string 
loop signal can be even larger than the amplitude for a 
string wake. However, both the amplitude and the size 
of the string-induced feature in the 21cm sky depends 
on G/i. Since the accretion pattern induced by a string 
loop is roughly spherical, the induced 21cm signal will 
be an ellipsoidal region of extra 21cm emission. Thus, 
even though the amplitude of the signal may be larger 
than that of a cosmic string wake, it will be harder to 
observationally distinguish from foreground noise. 



III. 21CM SIGNAL OF A COSMIC STRING 
LOOP 



A. Accretion of Matter onto a Cosmic String Loop 

We first study the accretion of matter onto a cosmic 
string loop. We consider a string loop which was pro- 
duced at time ti > t eq (and corresponding redshift zi). 
As a first approximation, we will replace the string loop 
of length I by a point source of mass m s = Ifi fixed at 
a point in space (which we take to be the center of the 
coordinates we use to study the accretion process). We 
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then consider the motion of mass shells surrounding the 
point mass shell in response to the gravitational shell, 
making use of the Zel'dovich approximation 27 1 . 

As just mentioned, we model the gravitational effects 
of the cosmic string loop in terms of a point mass with 
associated energy density 



P string = m s S(x) . 



(5) 



We will assume radial symmetry. In terms of the radial 
coordinate h the above energy density source becomes 



1 



Pstring Tils 



Anh 2 



5(h). 



(6) 



Let us now consider a mass shell surrounding the string 
loop. The physical distance of the shell from the mass 
point is 



h(q,t) = a(t)[q - ^(q,t)] , 



(7) 



where q is the comoving scale corresponding to the shell 
under consideration and ip(q, t) is the comoving displace- 
ment of this shell acquired in response to the gravita- 
tional force. The cosmological scale factor a(t) is nor- 
malized to be one today: a(io) = 1, where io stands for 
the present time. 

In the Zel'dovich approximation the dynamics of this 
mass shell is described by Newtonian gravity, i.e. by 



h = 



(8) 



where <E> is the Newtonian gravitational potential which 
in turn is determined by the Poisson equation 



V 2 = 4irG(p bg + Pstring) , 



(9) 



where pb g is the background energy density. 

Inserting the ansatz (J7J into the basic equations (|SJ) 
and © of the Zel'dovich approximation and linearizing 
in ip yields the following equation of motion for the co- 
moving displacement tp(q,t) 



eGtn 



qH 2 



(10) 



with initial conditions that both -0 and "0 vanish at the 
time U when the string loop is formed. This equation 
can be solved by means of a Born approximation. Since 
this is a standard calculation (see e.g. [3, H3|) we onr y 
quote the final result: 



1> = -. 



3 m q Gt 



3, t 



2 , t 



5 Ei 5 ti 



(11) 



The Zel'dovich approximation is good until the shells 
"turn around" , i.e. until h(q, t) = 0. This condition 
determines the comoving scale q n i(t) which is turning 
around at the time t. A simple calculation yields 

q nl = {\m s Gtlf'\l-fl^ 9o( £ )2/9 ' (12) 



where we have defined a quantity qo which we will use in 
the following. Note that this growth as a power of t 2 ^ 9 
agrees with the result of linear perturbation theory that 
the accreted mass (which is proportional to q^) scales as 
i 2 / 3 . 

After turnaround, the shell will collapse and virialize at 
a physical radius R max which is one quarter of the radius 
which the shell would have at the time of turnaround in 
the absence of gravitational accretion. Thus, the physical 
radius of the nonlinear structure seeded by the string loop 
at time t is 

R max (t) = \a(t)qni{t) = ^W" 1 ®^) 173 ( 13 ) 
and the total mass which is gravitationally bound is 



47T 

M(t) = —q nl (t) 3 p 0l 



(14) 



where po is the background energy density at the present 
time to ■ Inserting (|12[) and making use of the Friedmann 
equation to substitute for p yields 



M(t) 



' t n 2/3 



(15) 



in agreement with the result expected from linear per- 
turbation theory. 

The mean overdensity (in both baryons and dark mat- 
ter) in the nonlinear structure induced by a string loop 



p av = 64p 



(16) 



since the radius has shrunk by a factor of 4 compared to 
what it would be for unperturbed Hubble flow. In fact, 
towards the center of the structure the relative overden- 
sity in dark matter is higher than at the edges, since the 
shell which virializes at a distance R < R ma x from the 
center has virialized earlier when the background density 
was larger. The way we can compute the cold dark mat- 



ter distribution is to fix R < R„ 



determine the value 



of q for the shell which collapses to radius R, the time 
when that shell has collapsed, and from those data the 
mass M(R) within that shell. The result is M (R) ~ R 3/4 
which leads to a density scaling 



PDM (R) ~ ir 9 / 4 . 



(17) 



The details of this calculation are presented in [29( . They 
are the same which led to the conclusion that the velocity 
rotation curve about a cold dark matter halo is not flat - 
it is more peaked at the center - and that thus a pure cos- 
mic sting model of structure formation (which is ruled out 
because such a model does not yield the observed oscilla- 
tions in the angular power spectrum of CMB anisotropics 
[3(il ]) would require hot dark matter [3l|. We are inter- 
ested, however, in the density of baryons. Baryons will 
shock heat and their density distribution is expected to 
be uniform, with a relative overdensity given by (1161) . 
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We close this subsection by returning to the initial as- 
sumption we made, namely taking the entire loop mass 
to be concentrated at the origin. A necessary condition 
for this approximation to be justified is that the loop ra- 
dius is smaller than the physical height of the outmost 
virializcd shell. Replacing the loop radius by the loop 
length I yields a slightly more conservative condition 



: (t) > I. 



(18) 



For a loop formed at time U whose induced nonlinear 
structure is being considered at time t ^> U this condition 
becomes 



(19) 



and is hence easier to satisfy for small loops than large 
ones. Based on the scaling solution, these are the more 
numerous ones. Inserting numbers, if we are interested 
in considering string loops formed at t eq (corresponding 
to a redshift z eq ~ 3 x 10 3 ) which have the longest time 
to undergo gravitational accretion, and redshift z ~ 30, 
then the condition (|19[) is satisfied for strings with tension 
GVt = 1CT 7 for all I < U. 



B. Gas Temperature in the Nonlinear Structure 

We assume that the temperature of the gas inside the 
nonlinear structure induced by the cosmic string loop is 
characterized by the molecular kinetic energy which the 
particles obtain during the infall process. We use the 
virialization prescription that the infalling shell ends up 
with a height which is half of the height at turn-around, 
which is also one quarter of the height the shell would 
have in the absence of any gravitational accretion. If we 
denote the velocity of particles just before the shock as 
v shock , then the gas temperature is given by 



3/2k B T = l/2mv 



shock ' 



(20) 



where ks is the Boltzmann constant which we will usu- 
ally set to one in the natural units we are using. Thus, 
we must now compute v s h oc k- 

Based on the shell dynamics studied in the previous 
subsection, the physical velocity of the shall labelled by 
q at late times t (times when we can neglect the decaying 
mode of i/j(t)) is 



• _ 3 m s Gtl r 4 1 ^ 1/3l 
2 [ 5(t ti) 2 / 3 J ' 

The shock time t s hock is given by 



(21) 



(22) 



where hm.arXs h 



max 10 ' L max 



l/2a(t n i)q. We therefore have 



h(q,t sho ck) = ja{tni)l- 



(23) 



Based on our expressions for h(q,t) it is straightfor- 
ward to compute the time t s when the shell which turns 
around at time t n i(q) hits the shock. The result is that 
to a good approximation 



,2/3 _ 2/3, 1 + V2 
s ~ nl ( V2 ' ' 

Hence, the velocity at the shock is given by 



(24) 



,4/9 



-1/9 



Mi) 2 / 3 

Making use of (|20"]) . this leads to a gas temperature of 



(25) 



T K {t) = (l^ )mHI f/-)-^( ms Gtl)^ 

8/9 

1 +-2/9 



(26) 



M1) 4 / 3 

or, multiplying out the factors of order unity, 

T K (t) = [QA]m H i{m s G) 2 'h- A/9 t- 2 ^ . (27) 

The gas temperature Tk needs to be compared with 
the temperature T 7 of the background CMB photons to 
determine whether the 21cm signal from the string loop 
is in emission or in absorption. As a first step, we must 
insert the string loop mass m s into (f27j) . We express the 
string loop mass as a function of the ratio of the string 
length I to the formation time U: 



m s G = jGyd, l , 

£7. 



(28) 



and thus 

T K (t) = 0.4 x mHI (Gtf/ 3 (i ) 2/3 z(t)^zr ^ . (29) 

The background photon temperature is given by 

T 7 (t) ~ 3 x 10~ 13 GcVz(t) , (30) 
and hence the ratio of temperatures is 



mm 



TJt) 3 x 10- 13 GeV 



(G^/ 3 (f) 2/3 z(t)- 2 /V/ 3 



3,lO^G,)r(i) 2/ \(t)-^z-^ 

£7' 



(31) 



(where (G/x) 6 indicates the value of G[i in units of 10 -6 ). 
Evaluating this result for string loops formed at Zj = 10 3 
and observed at z 2 / 3 = 10 we obtain 



r 7 (t) 



3 x 10 6 a 2 / 3 (G M ) 2/3 



(32) 



6 



for the largest loops present at time ti and 



T K (t) 
T 7 (t) 



3 x 10 2 7 2 / 3 (G/x)^ /3 



(33) 



for loops which dominate the distribution at the time ti 
which have radius R = 'yG/iU at that time 5 

The main message from this subsection is that the gas 
temperature in the nonlinear structure formed by a string 
loop is typically larger than that of the CMB photons. 
Hence, the 21cm signal is an emission signal. This is a 
difference compared to the string wake case where 
for interesting values of Gfi the 21cm signal will be in 
absorption. The difference is due to the fact that string 
loops exert a stronger gravitational attractive force than 
that induced by string wakes. 



C. Brightness Temperature 

Consider the CMB radiation at frequency v in direc- 
tion of a gas cloud. The brightness temperature in this 
direction due to 21cm emission for a frequency v is de- 
noted by T b (y) and is given by (see [25j ) 



T b {y) =T s {l-e- T »)+T 1 (v)e 



(34) 



where T 7 is the background CMB temperature, t v is the 
optical depth and T$ the spin temperature of the gas. 
The hrst term describes the extra emission as a conse- 
quence of the hot gas, the second term represents the 
absorption of the background CMB radiation by the gas 
cloud. 

Because UV scattering is insignificant in our case, the 
spin temperature Ts depends only on the gas tempera- 
ture Tk of the wake and the collision coefficient x c : 



1 



1 



•T K 



(35) 



The collision coefficient x c describes the rate at which 
hydrogen atoms and electrons are scattered. We will use 
values of x c for which are of observational interest. The 
spin temperature T$ itself describes the relative number 
density 



n 



3exp(-T*/T s ) 



(36) 



of atoms in the two hyperfine states that when excited 
produce the 21cm radiation. The quantity T» = 0.0Q8K 



5 Loops with radius R < ati at time ti were formed earlier and so 
started accreting earlier if tj > t eq . For such loops (133 I t is not 
the correct formula to use. Instead, one should use the analog 
of l |32|l for the earlier initial time. However, since we are taking 
Zi ~ z eq . and small loops present at that time only start the 
gravitational accretion process at the time of equal matter and 
radiation, then 11331 is the equation which is applicable. 



is the temperature that corresponds to the energy split- 
ting between these two states. Furthermore, n\ and no 
are the individual number densities of atoms in the hyper- 
fine states. As before, t„ is the optical depth which can 
be determined by computing the absorption coefficient 
of the light way along the gas cloud. The frequency v 
is the blue shifted frequency at the position of the cloud 
corresponding to the observed frequency vq. The term 
proportional to T 7 is due to absorption and stimulated 
emission. 

As we have mentioned, what is of interest to us is the 
comparison of the temperature observed today in direc- 
tion of the hydrogen gas cloud of the defect as compared 
to the background temperature. This is given by 



ST b (u) 



T b {y)—T^{y) T s - T 7 



1 + z 



1 



(37) 



where we have assumed in the last step that the optical 
depth is smaller than 1 and that one thus Taylor expand 
the exponential factor to linear order in t„ 6 . The factor 
1 + z comes from the redshifting of temperatures between 
the time of emission and the present time. We can now 
express T 7 in terms of Tg to obtain 



6T b (y) ~ T s 



1 + X c 1 + Z 



(1 - TJT K ). 



(38) 



The optical depth for a general cloud of hydrogen is 
(in natural units) [25j 



3Ai N HI v w 



4^ 2 



4 T s 



(39) 



where A\q is the spontaneous emission coefficient of the 
21cm transition, Njji is the column density of hydro- 
gen within the gas cloud which the photons reaching us 
pass through, and <p(v) is the line profile. The hydrogen 
column density is given by the number density uhi of 
hydrogen atoms in the cloud multiplied by the physical 
diameter 2R of the cloud, where R is the height computed 
in the first subsection 7 : 



N. 



HI 



string 



R. 



(40) 



The line profile describes the broadening of the emission 
lines as a consequence of the spatial extent of the gas 
cloud and the resulting redshift difference 8v in the fre- 
quency 



6p 

— = 2HR, 
v 



(41) 



where H is the Hubble expansion constant. Because <p(v) 
is normalized to unity we find that 



(42) 



6 However in practice this result is skewed by the intcrgalactic 
medium which can also be considered a hydrogen gas cloud. 

7 Note that rifji is the hydrogen number density at the time t. 
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for v e [1/10 — vvq + Su/2] and 0(f) = otherwise. 

This frequency shift plays an integral role in the observ- 
ability of the 21cm signal since it determines the width 
of the string signal in the redshift direction. 

Note that the dependence on G\x cancels out between 
the column height and the line profile, as it did in the 
case of cosmic string wakes Q. Thus, we obtain 



.4 



10 



1 



riHi 



(43) 



where we have used H(z) = HoQm 2 {l + z) 3 ^ 2 in which 
Q m is the ratio between the matter energy density and 
the critical density (the density for a spatially flat uni- 
verse). We can insert this into the expression for 5T b (v) 
to obtain the signal strength 



ST b {u) 



x c 



1 + x, 



■(1 



T~ „ 3j4io 1 n 



string 
HI 



+ 



-5/2 



(44) 



We can remove the dependence on the hydrogen number 
density nui by introducing the background density n bg 
and noting that the overdensity of hydrogen n H r I m9 is 64 
times that of the background. This yields 



ST b (u) 



1 



■(1- 



T K '\Qv 2 



-4, 



string 



1/2 



(1+*) 



1/2 



(45) 

where we have rescaled the background number density 
to the present time. 

In the expression (|45j) for the brightness temperature, 
the string tension G\x enters only in two minor ways. 
Firstly, it determines the gas temperature Tk, but since 
Tjc 3> T-y for the parameter values we are interested in, 
this dependence is negligible. The second place where 
G[i enters is in the collision coefficient x c which depends 
on the gas temperature and hence on G/i. However, for 
our parameter values x c 3> 1 and hence this second de- 
pendence is negligible as well. 

Inserting the values Hq = 73kms~ 1 AIpc~ 1 , Am = 
2.85 x lO-^s" 1 , v w = U20MHz, O m = 0.26, T* = 
0.068A', and n bg = 1.9 x 10~ 7 cra~ 3 (l + z) 3 we obtain 



ST b (u) 



[1.11 



x c 



1 



1/2 



where the collision coefficient x c is given by 



,HH 

no 



AmT, 



10^7 



(46) 



(47) 



and Kjxf 1 is the de-excitation cross section whose value 
is 2.7 x 10 _9 m 3 s _1 at high temperatures corresponding 
to G/j, - 0.3 x 10~ 6 at a redshift of (1 + z) ~ 30 and 
(1+zi) ~ 10 3 . These values lead to the collision coefficent 
of x c » 1. We have also made use of the fact that the 
overdensity of the hydrogen gas is 64 times that of the 
background. The ensuing signal of the cosmic string loop 
is then 



&T h (y) ~GK, 



(48) 



which is an emission signal in the 21cm map. This tem- 
perature is even larger than that obtained for cosmic 
string wakes, the reason being that the overdensity in- 
side of a string loop-induced structure is larger than that 
in a string wake by a factor of 16. 



D. Geometry of the Signal 

Projected onto the two angular directions, the 21cm 
signal of a cosmic string loop looks like a filled circlular 
region of extra 21cm emission. The angular scale 9 of 
this region is determined by the ratio of the comoving 
distance £ com corresponding to R m ax and the distance of 
the observer's past light cone. We can express this as 

5SF-*f- (49) 

The comoving length is related to the physical distance 
by a(t): 

^corn & ^phys (1 ~\~ z)Rrnax j (^^) 

and thus the angular size of the defect is determined by 
the maximum radius of the nonlinear region which has 
collapsed onto the cosmic string loop. Thus, 

9 (l + z)R max {z) 



180° 



to 



(51) 



Starting from (fTB")) and inserting the mass of the string 
loop in terms of the loop length l s we obtain 



R-n 



^) = \{\f'\G,f' 3 ( l ff\z-^z-^ (52) 
The angular size therefore becomes 

9 ~ 180° X I(^)1/3 (Gm) 1/3(!£)1/3 z 1/6 z -1/3 . (53) 

4 5 t j 

Inserting l s = erf, (valid for the largest loops present at 
time tj), Zi = 10 3 and z — 30 we obtain 

9-2° xlQ- 1 (Gn)l /3 a 1/3 . (54) 

Smaller loops at ti give a correspondingly smaller angular 
scale. 

The redshift extent of the blob of extra 21cm emission 
coming from the nonlinear structure seeded by a string 
loop was derived previously. It is 



Sv 

— — 2HR Tl 
v 



(55) 



Inserting (|52[) and using the same parameters as above 
we obtain 



Sv 



2 x lQ- 6 (Gn) & ' a 



(56) 



which is a larger value than the thickness for the 21cm 
signal of a cosmic string wake. The reason for the larger 
value is that in the case of a string wake (planar accre- 
tion) the height of the nonlinear structure scales linearly 
in G/j, whereas for a string loop (spherical accretion) it 
scales as (G/^) 1 / 3 . 
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IV. DISCUSSION 

In this paper we have discussed the 21cm signature of 
cosmic strings loops. We found that for values of the 
string tension close to the current upper bounds, the sig- 
nature is an ellipsoidal region of extra 21cm emission. 
The brightness temperature of this signal is several de- 
grees K for string loops produced at a redshift of = 10 3 
and observed at a redshift of z = 30. The brightness tem- 
perature is to a first approximation independent of the 
string tension, as in the case of the signature of a cosmic 
string wake. 

Projected onto the celestial sphere, the angular extent 
of the region of extra 21cm emission is (for the same 
redshifts Zi and z mentioned in the previous paragraph) 
0.2° for a value of G[i = 1CU 6 and for a string loop with 
length l s = ti. This angle scale as (G^) 1 / 3 and also as 

1/3 

l s ' . In redshift direction the relative thickness of the 
region of extra 21m emission is — ~ 2 x 10~ 3 (for the 
same values of G/i and a as above), and this thickness 
has the same scaling in both G/i and l s as the angular 
scale. 

In terms of amplitude and both angular and redshift 
extent the cosmic string signal should be easily detectable 
by the SKA experiment, and possibly even with LOFAR. 
However, the signal might be hard to disentangle from 
noise and foregrounds, in contrast to the signal of a string 
wake which forms a very specific geometrical pattern. 

A scaling distribution of string loops will lead to a 
superposition of many ellipsoidal regions of extra 21cm 
emission. However, since 21cm redshift surveys will pro- 
duce three dimensional maps, the signatures of individ- 
ual string loops will be spaced further apart as the cor- 
responding signals in CMB anisotropy maps would be. 
Hence, we do not expect the string loop signal to become 
Gaussian as a consequence of the central limit theorem. 

Since the string loop distribution can be viewed as 
statistically independent from the distribution of string 
wakes, the 21cm signature of a scaling string model will 
be the linear superposition of what we found in the cur- 



rent paper with what was found in Q . It would be inter- 
esting to compute the total angular power spectrum of 
the 21cm signal, following what was done in (26|. How- 
ever, in computing a two dimensional correlation func- 
tion, and in doing this for a projected map, a lot of 
discriminatory power is lost. An improvement could be 
obtained by computing a three-dimensional power spec- 
trum. Even more promising, however, would be to find 
position space algorithms to search for string signatures. 

Our analysis is applicable with very few changes to 
a second type of topological defect, namely to global 
monopoles, pointlike defects arising in a theory in which 
a global symmetry is broken in a way that the second ho- 
motopy group of the vacuum manifold is non-trivial [32| . 
Note, however, that theories with broken global symme- 
tries have potential problems and are hence not consid- 
ered as interesting as those with broken local symmetries 
[33^ . Global monopoles also give rise to a spherical ac- 
cretion pattern and to nonlinear density perturbations 
at early times. The density distribution of a monopole 
is different from that of a cosmic string loop since field 
gradient energy extends arbitrarily far from the center 
of the monopole 8 . As a consequence, the dark mat- 
ter distribution has a different radial profile compared to 
that of the structure induced by a cosmic string loop. 
However, the baryon distribution inside the region which 
has undergone shock heating will be identical for string 
loops and global monopoles, and hence the 21cm signal 
of a single global monopole will not be distinguishable 
from that of a cosmic string loop (see [29| for details). 
On the other hand, the scaling distribution of monopoles 
will differ from that of string loops, and hence the overall 
21cm maps can in principle be distinguished. 



Acknowledgments 

This work is supported in part by a NSERC Discovery 
Grant and by funds from the CRC Program. 



[1] See |http:// www.lofar.org/ 1 . 

[2] C. L. Carilli and S. Rawlings, New Astronomy Review [7] 

48, 979 (2004). 
[3] R. H. Brandenberger, R. J. Danos, O. F. Hernandez and 

G. P. Holder, "The 21 cm Signature of Cosmic String [8] 

Wakes," JCAP 1012, 028 (2010) [arXiv: 1006.25141 [astro- 

ph.CO]]. 

[4] A. Vilenkin and E.P.S. Shellard, Cosmic Strings and [9] 
other Topological Defects (Cambridge Univ. Press, Cam- 
bridge, 1994). 

[5] M. B. Hindmarsh and T. W. B. Kibble, "Cos- 
mic strings," R ept. Prog. Phys. 58, 477 (1995) 
[arXiv:hep-ph/9411342l . 

[6] R. H. Brandenberger, "Topological defects and struc- [10] 
ture formation," Int. J. Mod. Phys. A 9, 2117 (1994) 



[arXiv:astro-ph/931004T] . 

A. H. Guth, "The Inflationary Universe: A Possible Solu- 
tion to the Horizon and Flatness Problems," Phys. Rev. 
D 23, 347 (1981). 

S. Sarangi and S. H. H. Tye, "Cosmic string production 
towards the end of brane inflation," Phys. Lett. B 536, 
185 (2002) [arXiv:hep-th/0204074] . 

R. Jeannerot, "A Supersymmetric SO(10) Model with 
Inflation and Cosmic Strings," Phys. Rev. D 53, 5426 
(1996) [arXiv:hep -ph/9509365] ; 

R. Jeannerot, J. Rocher and M. Sakellariadou, "How 
generic is cosmic string formation in SUSY GUTs," Phys. 
Rev. D 68, 103514 (2003) [arXiv:hep-ph/0308134l . 
E. Witten, "Cosmic Superstrings," Phys. Lett. B 153, 
243 (1985). 



9 



[11] R. H. Brandenberger and C. Vafa, "Superstrings In The 
Early Universe," Nucl. Phys. B 316, 391 (1989).; 
A. Nayeri, R. H. Brandenberger and C. Vafa, "Pro- 
ducing a scale-invariant spectrum of perturbations in a 
Hagedorn phase of string cosmology," Phys. Rev. Lett. 



97, 021302 (2006) arXiv:hep-th/0511140 



[13 



R. H. Brandenberger, A. Nayeri, S. P. Patil and C. Vafa, 
"String gas cosmology and structure formation," Int. J. 
Mod. Phys. A 22, 3621 (2007) [hep-th/0608T2l"] ; 
R. H. Brandenberger, "String; Gas Cosmology," 
larXiv:0808.0"746l [hep-th] . 
[12] C. Dvorkin, M. Wyman and W. Hu, "Cosmic String con- 
straints from WMAP and the South Pole Telescope," 
Phys. Rev. D 84, 123519 (2011) [arXiv: 1109.49471 [astro- 
ph.CO]]. 

L. Pogosian, S. H. H. Tye, I. Wasserman and M. Wyman, 
"Observational constraints on cosmic string pro- 
duction during brane inflation," Phys. Rev. D 68, 
023506 (2003) [Erratu m-ibid. D 73, 089904 (2006)] 
[arXiv:hep-th /0304188| ; 

M. Wyman, L. Pogosian and I. Wasserman, "Bounds on 
cosmic strings from WMAP and SDSS," Phys. Rev. D 
72, 023513 (2005) [Errat um-ibid. D 73, 089905 (2006)] 
[arXiv:astro-ph/0503364| ; 

A. A. Fraisse, "Limits on Defects Formation 
and Hybrid Inflationary Models with Three- Year 
WMAP Observations," JCAP 0703, 008 (2007) 
[arXiv:astro-ph/0603589] ; 

U. Seljak, A. Slosar and P. McDonald, "Cosmological 
parameters from combining the Lyman-alpha forest with 
CMB, galaxy clustering and SN constraints," JCAP 
0610, 014 (2006) [arXiviastro ph/0604335] ; 
R. A. Battye, B. Garbrecht and A. Moss, "Constraints 
on supersymmetric models of hybrid inflation," JCAP 
0609, 007 (2006) [arXiv:astro-ph/0607339l ; 
R. A. Battye, B. Garbrecht, A. Moss and H. Stoica, 
"Constraints on Brane Inflation and Cosmic Strings," 
JCAP 0801, 020 (2008) |arXiv:0710.154"T1 [astro-ph]]; 
N. Bevis, M. Hindmarsh, M. Kunz and J. Urrestilla, 
"CMB power spectrum contribution from cosmic 
strings using field-evolution simulations of the Abelian 
Higgs model," Phys. Rev. D 75, 065015 (2007) 



[arXiv:astro^pIi/ 0605018 ; 

N. Bevis, M. Hindmarsh, M. Kunz and J. Urrestilla, 
"Fitting CMB data with cosmic strings and inflation," 
Phys. Rev. Lett. 100, 021301 (2008) [astro-ph/0702223| 
[ASTRO-PH]]; 

R. Battye and A. Moss, "Updated constraints on the 
cosmic string ten sion," Phys. Rev. D 82, 023521 (2010) 
[arXiv: 1005.04791 [astro-ph.CO]]; 

J. Urrestilla, N. Bevis, M. Hindmarsh and M. Kunz, 
"Cosmic string parameter constraints and model analysis 
using small scale Cosmic Microwave Background data," 
JCAP 1112, 021 (2011) [arXiv:1108.2730l [astro-ph.CO]]. 

[14] S. Amsel, J. Berger and R. H. Brandenberger, "Detecting 
Cosmic Strings in the CMB with the Canny Algorithm," 
JCAP 0804, 015 (2008) |arXiv:0709.0982l [astro-ph]]. 

[15] A. Stewart and R. Brandenberger, "Edge Detection, Cos- 
mic Strings and the South Po le Telescope," JCAP 0902, 
009 (2009) [arXiv:0809.0"865l [astro-ph]]; 
R. J. Danos and R. H. Brandenberger, "Canny Al- 
gorithm, Cosmic Strings and the Cosmic Microwave 
Background," I nt. J. Mod. Phys. D 19, 183 (2010) 
[arXiv:0811. 20041 [astro-ph]]. 



[16] Y. B. Zeldovich, "Cosmological fluctuations produced 
near a singularity," Mon. Not. Roy. Astron. Soc. 192, 
663 (1980); 

A. Vilenkin, "Cosmological Density Fluctuations Pro- 
duced By Vacuum Strings," Phys. Rev. Lett. 46, 1169 
(1981) [Erratum-ibid. 46, 1496 (1981)]; 
T. W. B. Kibble, "Phase Transitions In The Early Uni- 
verse," Acta Phys. Polon. B 13, 723 (1982); 
T. W. B. Kibble, "Some Implications Of A Cosmological 
Phase Transition," Phys. Rept. 67, 183 (1980); 
N. Turok and R. H. Brandenberger, "Cosmic Strings And 
The Formation Of Galaxies And Clusters Of Galaxies," 
Phys. Rev. D 33, 2175 (1986); 

H. Sato, "Galaxy Formation by Cosmic Strings," Prog. 
Theor. Phys. 75, 1342 (1986); 

A. Stebbins, "Cosmic Strings and Cold Matter", Ap. J. 
(Lett.) 303, L21 (1986). 

[17] A. Albrecht and N. Turok, "Evolution Of Cosmic 
Strings," Phys. Rev. Lett. 54, 1868 (1985); 
D. P. Bennett and F. R. Bouchet, "Evidence For A Scal- 
ing Solution In Cosmic String Evolution," Phys. Rev. 
Lett. 60, 257 (1988); 

B. Allen and E. P. S. Shellard, "Cosmic String Evolu- 
tion: A Numerical Simulation," Phys. Rev. Lett. 64, 119 
(1990); 

C. Ringeval, M. Sakellariadou and F. Bouchet, "Cosmo- 
logical evolution of cosmic string loops," JCAP 0702, 
023 (2007) [arXiv:astro ph/051 1646] ; 

V. Vanchurin, K. D. Olum and A. Vilenkin, "Scaling of 

cosmic string loops," Phys. Rev. D 74, 063527 (2006) 

[arXiv:gr-qc/0511159] . 
[18] J. Polchinski and J. V. Rocha, "Analytic Study of Small 

Scale Structure on Cosmic Strings," Phys. Rev. D 74, 

083504 (2006) [arXiv:hep-ph/0606205l . 
[19] R. H. Brandenberger and N. Turok, "Fluctuations From 

Cosmic Strings And The Microwave Background," Phys. 

Rev. D 33, 2182 (1986). 
[20] A. Vilenkin, "Gravitational Field Of Vacuum Domain 

Walls And Strings," Phys. Rev. D 23, 852 (1981). 
[21] N. Kaiser and A. Stebbins, "Microwave Anisotropy Due 

To Cosmic Strings," Nature 310, 391 (1984). 
[22] J. Silk and A. Vilenkin, "Cosmic Strings And Galaxy 

Formation," Phys. Rev. Lett. 53, 1700 (1984); 

M. Rees, "Baryon concentrations in string wakes at z > 

200: implications for galaxy formation and large-scale 

structure," Mon. Not. R. astr. Soc. 222, 27p (1986); 

T. Vachaspati, "Cosmic Strings and the Large-Scale 

Structure of the Universe," Phys. Rev. Lett. 57, 1655 

(1986); 

A. Stebbins, S. Veeraraghavan, R. H. Brandenberger, 
J. Silk and N. Turok, "Cosmic String Wakes," Astrophys. 
J. 322, 1 (1987). 
[23] R. J. Danos, R. H. Brandenberger and G. Holder, "A 
Signature of Cosmic Strings Wakes in the CMB Polariza- 
tion," Phys. Rev. D 82, 023513 (2010) |arXiv: 1003.09051 
[astro-ph.CO]]. 

[24] R. Khatri and B. D. Wandelt, "Cosmic (super)string con- 
straints from 21 cm radiation," Phys. Rev. Lett. 100, 
091302 (2008) |arXiv:0801.4406l [astro-ph]]: 
A. Berndsen, L. Pogosian and M. Wyman, "Correlations 
between 21 cm Radiation and the CMB from Active 
Sources," larXiv: 1003.22141 [astro-ph.CO]. 

[25] S. Furlanetto, S. P. Oh and F. Briggs, "Cosmology 
at Low Frequencies: The 21 cm Transition and the 



10 



High-Reds hift Universe," Phys. Rept. 433, 181 (2006) 
|arXiv:astro-ph/0608032| . 

[26] O. F. Hernandez, Y. Wang, R. Brandenberger and 
J. Fong, "Angular 21 cm Power Spectrum of a Scaling 
Distrib ution of Cosmic S tring Wakes," JCAP 1108, 014 
(2011) |arXiv:1104.3337l [astro-ph.CO]]. 

[27] Y. B. Zeldovich, "Gravitational instability: An Approx- 
imate theory for large density perturbations," Astron. 
Astrophys. 5, 84 (1970). 

[28] R. H. Brandenberger, L. Perivolaropoulos and A. Steb- 
bins, "Cosmic Strings, Hot Dark Matter and the Large- 
Scale Structure of the Universe," Int. J. Mod. Phys. A 5, 
1633 (1990); 

L. Perivolaropoulos, R. H. Brandenberger and A. Steb- 
bins, "Dissipationless Clustering Of Neutrinos In Cosmic 
String Induced Wakes," Phys. Rev. D 41, 1764 (1990). 

[29] M. Pagano, MSc thesis, McGill University, 2012. 

[30] J. Magueijo, A. Albrecht, D. Coulson and P. Ferreira, 
"Doppler peaks from active perturbations," Phys. Rev. 
Lett. 76, 2617 (1996) |arXiv:astro-ph/9511042] ; 
U. L. Pen, U. Seljak and N. Turok, "Power 



spectra in global defect theories of cosmic struc- 
ture formation," Phys. Rev. Lett. 79, 1611 (1997) 
[arXiv:astro-ph/9704165l . 
[31] R. H. Brandenberger, N. Kaiser and N. Turok, "Dissipa- 
tionless Clustering Of Neutrinos Around A Cosmic String 
Loop," Phys. Rev. D 36, 2242 (1987); 
R. H. Brandenberger, N. Kaiser, D. N. Schramm and 
N. Turok, "Galaxy and Structure Formation with Hot 
Dark Matter and Cosmic Strings," Phys. Rev. Lett. 59, 
2371 (1987). 

[32] D. P. Bennett and S. H. Rhie, "Cosmological evolution of 

global monopoles and the origin of large scale structure," 

Phys. Rev. Lett. 65, 1709 (1990). 
[33] M. Kamionkowski and J. March-Russell, "Are tex- 

tures natural?," Phys. Rev. Lett. 69, 1485 (1992) 

[hep-th/9201063j ; 

R. Holman, S. D. H. Hsu, E. W. Kolb, R. Watkins and 
L. M. Widrow, "Cosmological texture is incompatible 
with Planck scale physics," Phys. Rev. Lett. 69, 1489 
(1992). 



